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Abstract 

We consider a set of M images, whose pixel intensities at a common point can be treated as 
the components of a M-dimensional vector. We are interested in the estimation of the mod- 
ulus of such a vector associated to a compact source. For instance, the detection/estimation 
of the polarized signal of compact sources immersed in a noisy background is relevant in 
some fields like Astrophysics. We develop two different techniques, one based on the 
Maximum Likelihood Estimator (MLE) applied to the modulus distribution, the modulus 
filter (ModF) and other based on prefiltering the components before fusion, the filtered 
fusion (FF), to deal with this problem. We present both methods in the general case of 
M images and apply them to the particular case of three images (linear plus circular po- 
larization). Numerical simulations have been performed to test these filters considering 
polarized compact sources immersed in stationary noise. The FF performs better than the 
ModF in terms of errors in the estimated amplitude and position of the source, especially in 
the low signal-to-noise case. We also compare both methods with the direct application of 
a matched filter (MF) on the polarization data. This last technique is clearly outperformed 
by the new methods. 

Keywords: Filters, image processing, matched filters, object detection, polarization, 
astronomy. 



1. Introduction 

The detection and estimation of the intensity of compact objects -i.e. signals with a 
compact support either in time or space domains- embedded in a background plus instru- 
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mental noise is a problem of interest in many different areas of science and engineering. 
A classic example is the detection of point-like extragalactic objects such as galaxies and 
galaxy clusters in sub-millimetric Astronomy. Regarding this particular field of interest, 
different techniques have proven useful in the literature. Some of the proposed techniques 
are frequentist, such as the standard matched filter 111], the matched multifilter [Q] or the 
recently developed matched matrix filters Q3|] . Other frequentist techniques include con- 
tinuous wavelets like the standard Mexican Hat [4] and other members of its family [|50 
and, more generally, filters based on the Neyman-Pearson approach using the distribution 
of maxima [@]. All these filters have been applied to real data of the Cosmic Microwave 
Background (CMB), like those obtained by the WMAP satellite CD and CMB simulated 
data [8D for the experiment on board the Planck satellite BSD, which has recently started 
operations. Besides, Bayesian methods have also been recently developed fio)] . Although 
we have chosen the particular case of sub-millimetric Astronomy as a means to illustrate 
the problem of compact source detection, the methods listed above are totally general and 
can be used to any analogous image processing problem. 

In most cases one is interested only in the intensity of the compact sources. In other 
cases, however, there are other properties of the signal that may be of interest. Such 
is the case, for example, of sources that emit electromagnetic radiation that is at least 
partially polarized. Polarization of light is conventionally described in terms of the Stokes 
parameters Q, U and V. Let us consider a monochromatic, plane electromagnetic wave 
propagating in the z-direction. The components of the wave's electric field vector at a 
given point in the space can be written as 

E x = a x (t) cos [aj t - 6 x (t)] 

E y = a y (t) cos [a) t-6y(t)y (1) 

If some correlation exists between the two components in the previous equation, then the 
wave is polarized. The Stokes parameters are defined as the time averages 

/ = (a 2 x ) + (a 2 y ), 

Q = (a 2 x )-{a 2 y ), 

U = (2a x a y cos [6 X - 6 y \), 

V = (2a x a y sm(e x -6 y )). (2) 

The parameter / gives the intensity of the radiation which is always positive, while the 
other three parameters define the polarization state of the wave and can have either sign. 
Q and U are the linear polarization parameters and V indicates the circular polarization of 
the wave. Unpolarized radiation is described by Q = U = V = 0. 
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While the total intensity of the wave is independent of the orientation of the x and y 
axes, the values of the other Stokes parameters are not invariant with respect to changes of 
the orientation of the receivers. On the other hand, the total polarization, defined as 



P = VG 2 + U 2 + V 2 , (3) 

is invariant with respect to the relative orientation of the receivers and the direction of the 
incoming light, and therefore it is a quantity with a clear physical meaning. For the case of 
purely linear polarization, V = and the previous expression reduces to P = (Q 2 + U 2 ) l/2 . 
Note that in order to get P from its components Q, U, V it is necessary to perform a 
non-linear operation. 

Although strictly speaking Q, U and V are the components of a tensor, the invariant 
combination © gives a quantity that can be seen as the modulus of a vector. In this paper 
we will introduce a methodology that can be applied to any problem in which a set of im- 
ages contain signals whose individual intensities can be considered as the components of 
a vector, but where the quantity of interest is the modulus of such a vector. For illustrative 
purposes, throughout the paper we will use as an example the case of light polarization, but 
the methods we will introduce are not limited to the example. Another possible application 
could be the determination of the modulus of a complex-valued signal. For example, in 



f| 1 If! the case of the estimation of the modulus of a complex- valued quantity whose com- 



ponents follow a Gaussian distribution was addressed. The techniques presented in their 



paper are related to our methods, but restricted to the two-dimensional case. In II 1211 . four 



methods: MLE estimator, median estimator, mean estimator and Wardle and Kronberg's 



estimator fl 1 311 are applied to the estimation of polarization (two-dimensional case). How- 
ever, these methods are applied in the cited paper to a single data and cannot be, except for 
the MLE, generalized to the detection of a signal with a given profile in a pixelized image, 
as considered in our paper. This is due to the fact that these methods lead to a system of 
possibly incompatible equations. 

Going back to the case of point-like extragalactic objects in sub-millimetric Astron- 
omy, the polarization of the sources plays and important role in the cosmological tests de- 
rived from CMB observations. Standard cosmological models predict that CMB radiation 
is linearly polarized. However, some cosmological models predict in addition a possible 



circular polarization of CMB radiation II 1411 . In order to better constrain the cosmological 
model with observations, it is crucial to determine the degree of polarization of not only 
the CMB radiation but also the other astrophysical sources whose signals are mixed with 



it. For an excellent review on CMB polarization, see Ill5il . We have treated the application 



to the detection of linearly polarized sources in CMB maps elsewhere lllql . but it is known 



that extragalactic radio sources can also indeed show circular polarization Ill7ll . Besides, 



circular polarization occurs in many other astrophysical areas, from Solar Physics fl 1 811 to 
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interstellar medium [19|], just to put a few examples. Therefore, in this paper we aim to 
address the general case of linear plus circular polarization of compact sources. 

In the general case we have three images Q, U and V. Different approaches can be used 
to deal with detection/estimation of point-like sources embedded in a noisy background. 
On the one hand, one can try to get the source polarization amplitude A directly on the 
P-map. In this approach, we will consider one filter, obtained through the MLE applied to 
the modulus distribution (ModF). On the other hand, we can operate with three matched 
filters, each one on Q, U and V followed by a quadratic fusion and square root (FF). We are 
trying to compare the performance of the two techniques for estimating the position and 
polarization amplitude of a compact source. Of course, in the case we have only the map 
of the modulus of a vector and the components are unknown, the FF cannot be applied 
but we can still use the ModF. Finally, we also apply a matched filter (MF) directly on 
the P-map and compare this simple method with the two new methods introduced in this 
paper. 

In Section[2]we will develop the methodology for the case of M images, because of the 
possible interesting applications to the general M-dimensional case and in particular to the 
3-dimensional case (polarization). In Section [3] we will show the results when applying 
these techniques to numerical simulations of images of Q, U and V that are relevant for 
the detection of compact polarized sources in Astrophysics. Finally, in Section |4] we give 
the main conclusions. 



2. Methodology 

2.1. The case of M images 

To develop our methodology, let us consider M images with intensities dj(x) at each 
point x of their common domain, j = 1, . . . ,M. Hidden in these images there is an un- 
known number of signals. In this work we consider signals with compact support ('com- 
pact sources' hereinafter), as for example galaxies in CMB polarization images. For sim- 
plicity, let us consider the case of a single compact source embedded in the images. We 
will assume a linear model 

dj(x) = Ajt(x) + tij(x). (4) 

In this equation, the source is characterized by amplitudes Aj in each image and by a spa- 
tial profile t(x) that is the same for the M images. This last condition is satisfied when 
the instrument resolution is much higher than the source scale . The source is immersed 
in noise nj(x) that is Gaussian and independently distributed with zero mean and disper- 
sion o-(x). Then the distribution of dj(x) is Gaussian with mean Ajt(x) and dispersion 
cr{x). We will consider that the noise is non- stationary but the dispersion is the same for 
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the different images at the same pixel. These conditions are typical in polarization im- 
ages. By construction, the total polarization map P(x) = (zZjdj(x)) l/2 includes a source 
characterized by a total polarization amplitude A = (£ ; Ap 1/2 , the modulus of the vector 
(A l ,A 2 ,....,A M ). 

2.1.1. ModF on the P-map 

If the noise is distributed normally and independently, the 1-pdf of dj at any point of 
each image is 

'-(dj(£)-AjT(Jt)) 2 



f(dj(x)\Aj) = 



1 



V2ncr(x) 



exp 



2(T 2 (x) 



(5) 



where t(x) is the known source profile at the corresponding point. Since the noise is 
independent between the different images, the M-pdf can be written as 



f(di(x), ...,d M (x)\A u ...,A M ) 



1 



(2n) M/2 o- M (x) 



exp 



(dj(X)-AjT(x)) 2 

2cr 2 (x) 



(6) 



Next, we will derive the distribution of P{x), the modulus of {d\{x), d 2 (x), ....,d M {x)). By 

changing to M-dimensional spherical coordinates, (P, 6>i, , 6 M -\), we write at any point 

x, (in order to simplify the notation we write P instead of P(x) ) 



jM-l 



f{p,e u ...,8 M -i\A) 



(2tt) 



exp 



P 2 +A 2 t 2 + PAt cos e M -\ 



2o- 2 



X 



(sin0 M _!) 



M-2 



(sin^ 2 ) M - 3 ...sin^ 2 



(7) 



This is the joint pdf of the modulus P and the corresponding angles, note that A = 
(zZj-Aj) 1 ^ 2 and we have multiplied by the Jacobian of the coordinate change. Finally, by 
integrating on the angles 8\, , 6m-u we find the pdf of the modulus P: 



f(P\A) 



pM/2 



exp [-{A 2 t 2 + P 2 )/2o- 2 ]Im_2 (at-^) , 



(8) 



-2( At )(M-2)/2 

where 7(m-2)/2 is the modified Bessel function of the corresponding order. When M=2 we 
obtain the Rice distribution rfiofr . In the case cr = 1 with a general M, distribution ([8]) is 
the non-central chi distribution with M degrees of freedom and A = At 112 U 12211 . This was 
expected since we have derived the pdf of the modulus of a vector whose components fol- 
low independent Gaussian distributions with mean Ajt. If there is no source, the previous 
formula defaults to 

/™= wr(1+M/2) exp[- f V2^]. (9) 
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where Y is the gamma function. In the case cr = 1, Q is the chi distribution. In particular, 
when M=2 we obtain the Rayleigh distribution 12311 . 

Now we assume that the M images are pixelized with the same pixel size; since the 
noise is independent pixel to pixel, the different values of P at each pixel, P h i = 1, . . . , N, 
with N the number of pixels, follow the two distributions 

MP M 1 

f(P u .., P N \0) = f] 2 M /Vfr(M/2 + 1} [-Ph2aj] (H ), (10) 



,M/2 



f(P 1 ,...,P N \A) = \ 



exp 



A 2 rj+Pj 
2o- 2 



crj(ATd 2 



PiTi 



(ID 



being H a and Hi the null (absence of source) and the alternative (presence of source) 
hypotheses, respectively, and r, the profile at the z th pixel. The log-likelihood is defined by 



~ a2 Z ^ " N(M ~ 2)12 logA + E log /(m ~ 2)/2 ( A ~?) 



/(A|P 1 ,...,P iV ) = log/(H 1 ) = 



(12) 



Where we have only written the terms wich depend on A. The MLE of the amplitude, A, 
can be obtained by maximizing the previous expression and is given by the solution of the 
equation 

yt = —■ (13) 



1^-2=^1 

: °t : i(MI2-\) 



(Ay,) 



(Ayd 



This equation can be interpreted as a non-linear filter operating on the data (ModF). This 
method allows us to obtain the MLE of A given a total polarization map P. 

2.1.2. Filtered Fusion ( FF) 

In this case we consider separately each one of the M images. We call dp the intensity 
in the j th image at the / th pixel and cr, the dispersion at each pixel. The pdf of the intensities 
in the j' th image is 



f(d fl , ...,d jN \Aj) = 

The MLE for Aj yields 
* _ LidjiTjlaj 

1 ~ 2/2 



1 



(2^) JV/ V 1 a N 



exp 



Z 



(dp - AjTi) 2 
2o-J 



(14) 



(15) 
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This is equivalent to the application of the linear matched filter MF, O, operating on each 



image dj, as given by Q24J] 



O(f) oc J^L. (16) 

Then, with the M filtered images A , we make the non-linear fusion A = ( £ • Ay 1 , in order 
to estimate the source polarization amplitude. We call this operation FF. Note that this 
method is clearly different from the ModF, now we calculate the MLE of each component 
Aj, given the values dji of that image and finally compute the modulus A of that estimated 
vector, whereas with the ModF we calculate the MLE of A given the polarization map P. 
For using the FF, we need to know all the images. On the other hand, for using the ModF 
we only need to know the total polarization map P. 

2.1.3. Matched Fiter (MF) 

This is a naive approach, assuming that we only know the P-map. We estimate A 
according to the expression 

« liiPiTi/aj 

A = v 2 2 < 17 > 

with Pj the polarization data at each pixel. This simple technique can be useful for com- 
parison with the more elaborated methods presented in the previous subsections. 

2.2. The case of three images 

Now, we shall assume that we have the same compact source in three images Q, U, 
V (this is the standard notation for Stokes parameters), characterized by amplitudes Aq, 
A v , A v and a profile t(x) immersed in noise nQ UV (x) that is Gaussian and independently 
distributed with zero mean and dispersion cr(x). In general, we will consider that the noise 
is non-stationary. We will assume a linear model for the three images 

Q, U, V(f) = A Q>[/iV T(f) + n QMV (^). (18) 

The P-map, P(x) = (Q 2 (x) + U 2 (x) + V 2 (x)) l/2 , is characterised by a source with amplitude 
A = (A 2 Q +A 2 U+ A 2 V )^ 2 . 
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2.2.1. ModF on the P-map 

In the case of three images, formulas (PT21) and (fl3l) can be written as 

1{A\P U ...,P N ) = log /(Hi) = -A 2 £ ^3 -NlogA + ^log sinhU^] 



(19) 



-i + T = yy i coth(A 3 ; ! ), yi=-T- (20) 
(T A ^ erf 



This equation can be interpreted as a non-linear filter operating on the data which is the 
(ModF) for this particular case. 

2.2.2. Filtered fusion (FF) 

In this case, we use the same MF operating on each image Q, U, V, as given by 
equation (fT5T ). Then, with the three filtered images Q M f, U M f, Vmf we make the non- 
linear fusion A = (Q 2 MF + U 2 MF + V 2 MF ) 111 pixel by pixel. 

2. 2. 3. Matched Filter (MF) 

We just apply (fT71) to the case of three images. 

3. Simulations and results 

As commented in the introduction, since the estimation of the intensity of polarized 
sources is of great interest in Astrophysics, we will use an example taken from CMB As- 
tronomy in order to illustrate the performance of the techniques introduced in the previous 
section. 

In order to compare and evaluate the performance of the two filters, we have simulated 
images of 16 x 16 pixels with a pixel angular sizeQ of 3 arcmin. We simulate the Q, U 
and V components of the polarization as follows: each component consists of Gaussian 
uncorrected noise, plus a polarized point source filtered with a Gaussian-shaped beam 
whose full width half maximum (FWHM) is 14 arcmin. This is a typical example of a 
CMB polarization experiment So the source polarization components can be written as 



s q,u,v = Ag^yexp 



'ly 2 



(21) 



'The angle of the sky subtended by a pixel of the detector, for a given telescope. 

2 This particular choice of the pixel and beam sizes corresponds to the specifications of the 70 GHz 
channel of the ESA's Planck satellite. 
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where y is the beam dispersion (angular size) and we assume in this formula that the source 
is centered at the origin. We consider stationary noise,with zero mean and r.m.s. deviation 
cr = 1 in some unit system^. We take values of Aq, A v and A v ranging from to 2.5 with 
a step of 0.5. The number of simulations is 100 for each combination of triplets of values 
of Aq, Au and A v . After carrying out the corresponding simulations for Q, U and V, we 
add them quadratically and take the square root to calculate P = yjQ 2 + U 2 + V 2 , the total 
polarization. 

We assume that we do not know the exact position of the source in the map and then 
we place it at random in the image. We have considered images of 16 x 16 pixels in order 
to do fast calculations. In order to avoid border effects, we simulate and filter 24 x 24 pixel 
patches and, after the filtering step, we retain only the 16x16 pixel central square. 

We use three different filters: the FF, as described in section I2~. 2 .21) . which consists in 
the application of the matched filter to the images in Q, U and V separately and then the 
calculation/construction of the P-map from the matched- filtered images, the ModF applied 
directly on P, derived from the MLE applied to the modulus distribution and presented 
with detail in section I2~.2. II and a simple MF applied on the P-map. We apply these filters 
to each simulation, centering the filters successively at each pixel, since we do not know 
the source position. We estimate the source amplitude A ModF for the ModF, in this case 
we calculate the value of A which maximises the log-likelihood, eq. (fT9~l) . For the FF we 

estimate separately and obtain A FF = ^Q 2 MF + U 2 MF + V^ F . With the MF, we obtain the 
estimator of A, (TT71) . Then, we have constructed three maps ( one for the ModF , one for 
the FF and another one for the MF) with the estimated values of A at each pixel. 

We show in Figure[T]four images corresponding to a polarized source with (Aq,Au,Av) = 
(1.5, 1.5, 1.5) embedded in noise. For the sake of a better visualisation, we show 48 x 48 
pixel images instead of the 16x16 sized images used in the simulations. We show the 
original image in P including noise and source, the image of the source, the image filtered 
with the ModF and finally, the image treated with the FF method. As we will comment 
with more detail below, the performance of the MF is much worse than that of the other 
filters and we will only show some results obtained with the MF, leaving the figures and 
tables for the comparison of the FF and the ModF. 

We compute the absolute maximum of each filtered map {A ModF , A FF and A M F) and 
keep this value as the estimated value of the polarization amplitude of the source and the 
position of the maximum as the position of the source. Note that for the more realistic case 
where more than one source can be present in the images, it is still possible to proceed as 



3 For this example, we use arbitrary intensity units, since the quantity of interest for our purposes is the 
signal to noise ratio of the sources. 
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20 40 20 40 



Figure 1: Top left: Image of a polarized source filtered with a Gaussian beam. The source polarization 
components are {Aq,Au,Av) = (1.5, 1.5, 1.5) and it is embedded in stationary noise <x = 1 . Top right: 
Image of the polarized source only. Bottom left: The first image after application of the ModF. Bottom 
right: The first image after the application of the FF. 

described by looking for local peaks in the filtered images. 

We also calculate the significance level of each detection. In order to do this, we carry 
out 1000 simulations with Aq = 0, A v = 0, A v = and we calculate the estimated value of 
the source polarization in this case for each filter. We consider the null hypothesis H (there 
is no polarized source) against the alternative hypothesis Hi (there is a polarized source). 
We set a significance level a = 0.05; this means that we reject the null hypothesis when 
a simulation has a estimated source polarization amplitude higher than that of 95% of the 
simulations without polarized source. The previous significance fixes a lower threshold 
(A* = 1.99 for the ModF and A* = 1.17 for the FF and A* = 3.81 for the MF) defining a 
region of acceptance in the space of polarization amplitudes. We define the power of the 
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test as 1 - 6, with 6 the probability of accepting the null hypothesis when it is false, i.e. the 
power is the proportion of simulations with polarized source with a estimated amplitude 
higher than that of A*. The higher the power the more efficient the filter is for detection. 

We calculate the estimated value of the polarization amplitude A and compare it with 
the real value A, obtaining the relative error of the estimation and its absolute value for 
each simulation.We also calculate the estimated position of the detected source and obtain 
the position error expressed in terms of the number of pixels. The average of all these 
quantities for 100 simulations are presented in Table Q] for all the cases. The rows in the 
Table are sorted in ascending order of A. We see in the table that the power for the FF 
is higher than for the ModF. The improvement is particularly high for A values equal or 
lower than 2.5. In this case the position and polarization errors are also lower for the FF. 
The FF can detect sources from A > 1.8 with power > 0.99 and average relative error 
(bias) < 0.05, average of its absolute value < 0.14 and average position error < 0.53. For 
A = 1.8, the ModF detects with power = 0.41, average relative error = 0.25, average of 
its absolute value = 0.25 and average position error = 1.39. For A = 1.8, the MF detects 
with power = 0.46, average relative error = 1.20, average of its absolute value = 1.20 
and average position error = 1.34. The errors of the polarization estimation for the MF 
are much higher than for the other methods, whereas the power and the position error are 
similar to those obtained with the ModF. 

The ModF and the FF perform in a similar way for signal-to noise ratio A > 3, in this 
case the power is 1 for both filters and the bias is < 0.02. The ModF is also slower than 
the FF, due to the maximization process involved in its application. However, note that the 
FF cannot be applied when we only have the image of the modulus i.e. we do not have 
information about the components, in this case the ModF can be a suitable filter, specially 
for signal-to-noise ratio > 3. 

In Figure [2} we have plotted the estimated polarization amplitude A against the real 
polarization amplitude A for the ModF and the FF. The average and 68% confidence in- 
tervals of 100 simulations are shown. It is clear that the FF performs better till A = 3. 
In Figure [Q we can also see the better performance of the FF. From a qualitative point 
of view, the ModF image shows more structure, with bright artifacts that could lead to 
spurious detections, while the FF image looks smoother. This is easy to understand if we 
think ModF is enhancing an inherently non-Gaussian noise, whereas FF is the (non linear) 
composition of three smoothed Gaussian noises. 

Finally, we have plotted in Figure [3] the position error in numbers of pixels for 1000 
simulations in the particular case Aq = 0.5, Ay = 1, Ay = 1. Only the simulations with 
a significance a = 0.05 have been taken into account and represented. It is clear that the 
position error is lower for the FF. 

In order to try to understand the better performance of the FF, we have carried out 
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1 2 3 4 5 

A 




Figure 2: Estimated source polarization amplitude A plotted against the real polarization amplitude A. The 
average and 68% confidence intervals of 100 simulations are plotted. Top : the ModF has been used. Bottom: 
the FF has been applied. The straight line A - A has been drawn for comparison. 

simulations of three-dimensional vectors whose components are gaussian-distributed with 
dispersion <x = 1 and values of the mean ranging from 0.5 to 2.5. We have considered 
the ModF, i.e we find the value of A which maximizes © with r = 1 and the filtered 
fusion, which amounts to averaging the components and calculating the modulus of this 
average vector; this last method is also considered in [1 lj] for the 2-dimensional case. The 
conclusion for these simple cases is that the FF performs better than the ModF, specially 
for the low signal-to-noise case, i.e. these simple examples confirm our conclusions: it is 
more precise to estimate the modulus estimating first the components and then calculating 
the square root of the quadratic sum than calculating the MLE estimator directly on the 
modulus. 



12 



15 



10 



O 



> 



-5 



-10 



o 



o 



* 



* 


ModF 


o 


FF 


* 




* 




* 




* 





o 



* 
o o 



o ® ® ® o 

iT*> ^T>. <g?> <Oi /CDv 
v©7 ^7 

*/<ps Ad /®j jj^. 
\£7 

* * * 

* 



* 



O 



* 



* * * 

(OX \U \U_ 



o 
* * 



* * 



* * 

O 



* 



-15 



-15 



-10 



-5 





X 



10 



15 



Figure 3: Position error in number of pixels for the case Aq — 0.5, Ajj - I, Ay - 1. 1000 simulations have 
been considered. 



4. Conclusions 

In this paper, we deal with the detection and estimation of the modulus of a vector, 
a problem of great interest in general and in particular in astrophysics when we consider 
the polarization of the cosmic microwave background (CMB), extragalactic sources, the 
interestellar medium or the Sun. The total polarization intensity P is defined as P = 
(Q 2 + U 2 + V 2 ) l/2 , where Q , U ( linear polarization) and V (circular polarization) are the 
Stokes parameters. We consider the case of images in Q, U and V consisting of a compact 
source with a profile t(x) immersed in Gaussian uncorrelated noise. We intend to detect 
the source and estimate its polarization amplitude by using three different methods, two 
new and a standard matched filter, a) a filter operating on the modulus P (ModF) and 
based on the maximization of the corresponding log-likelihood and b) a filtered fusion 
(FF) procedure, i.e. the application of the matched filter (MF) on the images of Q, U and 
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V and the combination of the corresponding estimates by making the non-linear fusion 
A = (Q 2 MF + U 2 MF + V 2 MF ) XI2 and c) a matched filter (MF) operating on the P-map. We 
present the three filters in Section[2]for the general case of the modulus of a M-vector, and 
for the case of three-dimensional vectors, deriving the corresponding expressions , eqs. 
([T3l and (1201) . for the estimation of the modulus of a vector with the ModF. 

Since we are interested in the detection of polarized signals in Astrophysics, we have 
only considered the three-dimensional vector case in our simulations. Note, however that 
the methods can be applied in the general case of combination of M images to obtain the 
image of the modulus of a vector. 

We have compared the performance of the filters when applied to simulated images 
consisting of the Q, U and V components of a Gaussian- shaped signal with different in- 
tensities plus Gaussian uncorrelated stationary noise. For each simulation and for the two 
methods, we have estimated the source amplitude and position. Besides, we have calcu- 
lated the detection power for a fixed significance a = 0.05 and the errors of the estimated 
amplitude and position. We find that the performance of the FF is the best for low signal- 
to-noise, the ModF performs like the FF for A > 3. The MF produces very high errors in 
the polarization estimation, making this filter unsuitable for the treated problem. 

We want to point out the good performance of the ModF, which could be an interesting 
alternative to the FF when we have an image of the modulus, but we do not know the 
components of the vector. 

Finally, it would be interesting to generalize our methods to the case of images in 
which the noise at different pixels is not uniform or, more generally, is correlated. This 
could be important for CMB polarization observations. We leave this problem for further 
work. 
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